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I. INTRODUCTION 



The solutions of the Dirac equation having the pseudospin, and spin symmetry have been 
strongly studied in the last years. The concept of the pseudospin symmetry [1, 2, 3] is an 
important theme in nuclear theory because of its features related to construct an effective 
shell-model coupling scheme [4], to study the structures of the deformed nuclei [5, 6]. The 
pseudopin symmetry occurs in nuclei when the magnitude of the scalar and vector potentials 
are nearly equal, but opposite sign, i.e., Vv{r) ~ — K(r) ,where the scalar potential is negative 
(attractive) and the vector potential is positive (repulsive) in the relativistic region. Further, 
the spin symmetry appears when the magnitude of the scalar and vector potentials are nearly 
equal, i.e., Vy{r) ~ Vs{r) [7-12]. The Dirac equation under the exact pseudospin and/or spin 
symmetry has been studied by with different type of potentials such as the Hulthen potential 
[13], the Morse potential [14], the Woods-Saxon potential [15, 16], and harmonic oscillator 
[17-21]. 

In the present work, we deal with the solutions of the Dirac equation if the exact pseu- 
dospin, and spin symmetry occur in the theory under the effect of a class of potentials, which 
have exponential form, i.e. the hypergeometric-type potential [22], the generalized Morse 
potential and the Poschl-Teller potential [23]. 

In order to obtain the energy eigenvalue equation and the corresponding wave functions, 
we apply a new approximation scheme [24], which is the parametric generalization of the 
Nikiforov-Uvarov (NU) method [25], by using an approximation to the centrifugal-like term. 
So, we obtain the energy spectra of the above potentials for the spin-orbit quantum number 
K — 0, or for any value for the case of pseudospin and spin symmetry, respectively. 

The organization of the present work is as follows. Firstly, we give briefly the equations 
for the Dirac spinors including the centrifugal term, and the spin-orbit quantum number k. 
Than we present the basics of the parametric generalization of the NU method. We find the 
bound states and the corresponding wave functions of the Dirac equation with the above 
potentials in the case of pseudospin and spin symmetry, respectively. Finally, we give our 
conclusions. 
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II. DIRAC EQUATION 



The Dirac equation for a particle with rest mass m in the absence of the scalar and vector 
potential is written {h — c— 1) [15] 



{a.p + I3[m + K(r)]}*(r) = [E - K(r))*(r) , (1) 
where p is the momentum operator, 5, and j3 are 4x4 Dirac matrices, i.e. a = 



^ a,: 



/3 = 



h 



V 





/ is the 2x2 unit matrix, and ai{i — 1, 2, 3) are the Pauh matrices. The 



V 



-I 

operator K is the spin-orbit matrix operator and written in terms of the orbital angular 
momentum operator L a,s K — —^(a.L + 1) ,which commute with the Dirac Hamiltonian. 
The Dirac spinors can be labelled by the quantum number set {n, k) ,where k is the eigenvalue 
of the spin-orbit operator, and written as 



*n«(r) = r "'^ h (2) 




where = [F„«(r)/r]y/^(^, 0) is the upper, and g^^ = [iG„„(r)/r]y/^(^, 0) is the lower 
component, and Yj^{9, 0) , and Yj^{6, (j)) are the spherical harmonics, respectively. The total 
angular momentum, the orbital angular momentum, and pseudo-orbital angular momentum 
can be written in terms of the spin-orbit quantum number /t = ±1, ±2, . . . , such as j = 
\k\ - 1/2 , i = \k + l/2\ - 1/2 , and e = \k - l/2\ - 1/2, respectively. 

Substituting Eq. (2) into Eq. (1), and eliminating ^^^(r), we obtain two uncoupled 
differential equations for the lower, and upper components of the Dirac equation 

f (P k(k — 1) 1 dE(r) / d ^^\^ ^ , . 



and 



f (P k(k + 1) 1 dA(r) / d k\i,, , ^ , . , ^ , . ^ , . 

{j^. - ^ + + 7)}^"M = M.WM.WF„M . (4) 

where MA(r) ^ m + E - A(r) , A(r) = K(r) - Vs{r) , Mj:{r) ^ m - E + E(r) , and 
S(r) = K(r) + K(r). 
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III. PARAMETRIC GENERALIZATION OF THE METHOD 



Let us give briefly the parametric generalization of the NU method. By using an ap- 
propriate coordinate transformation, the Schrodinger equation can be transformed into the 
following form 



(5) 



dz"^ ' ^ ' ^ ' dz 

where (t{z) , and a{z) are polynomials, at most, second degree, and f{z) is a first-degree 
polynomial. By writing the general solution as ^(-z) = ip{z)(fi{z), we obtain a hypergeometric 
type equation [24] 



dz^ a{z) dz 
where ip{z) and (pn{^) are defined as [24] 



a{z) 



(6) 



1 dip{z) Tr{z) 



ip{z) dz a{z) ' 



(7) 



Vn{z)^^^K{z)p{z)], 



(8) 



p(z) dz" 

where a„ is a normalization constant, and p{z) is the weight function satisfying the following 
equation [24] 



dz p{z) dz 

The function 7r{z) , and the parameter A in the above equation are defined as 



(9) 



n{z) = - [a'{z) - f{z)] ± 
X = k + n'{z) . 



-[a'iz)-f(z)r-a(z) + ka(z) 



1/2 



(10) 
(11) 



In the NU method, the square root in Eq. (10) must be the square of the polynomial, so 
the parameter k can be determined. Thus, a new eigenvalue equation becomes 



X^X„^-nT'{z)--{n'-n)a"{z) 



(12) 



where prime denotes the derivative and the derivative of the function t{z) = f{z) + 27r{z) 
should be negative. 

Now, in order to clarify the parametric generalization of the NU method [25] , let us take 
the following Schrodinger-like equation written for any potential 



2 .^d^^iz) 
z {1- asz) 



+z{l-a3z){ai-a2z)^ + [-Ciz''+C2Z-C3]'^{z)^0. (13) 



Comparing Eq. (13) with Eq. (5), we obtain 



f{z) = «! - a2Z ; a{z) = z{l - a^z) ; a{z) = -^iz'^ + - 6 
Substituting these into Eq. (10), we obtain 



(14) 



7r{z) = a4 + a^z ± 



with the following parameters 



(cKg — kasjz"^ + {aj + k)z + ctg 



1/2 



(15) 



(16) 



q;4 = 1(1 - ai) , "5 = 1 («2 - 203) , 
aQ — a\-\-^i, — 2a4a5 — ^2 , 

as = + 6 ■ 

We obtain the parameter k from the condition that the function under the square root 
should be the square of a polynomial 



ki,2 = -(a? + Sttattg) ± 2^/a^, 
where ag = a^a-j + a^as + . The function t:{z) becomes 



(17) 



'k{z) = q;4 + ar,z - [{^fa^ + a^^fa^ )z - y^] 



(18) 
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for the /c-value k — — (ccy + 2a3as) — 2^/0309 . We also have from t{z) — f{z) + 2Tr{z) , 

r{z) = cci + 2a4 — (0:2 - 20:5)2; - 2 [{^/a^ + as^/a^)z — ^/a^] . (19) 
Thus, we impose the following condition to fix the /c-value 

t'{z) = -{a2-2a5) -2{^/a^ + a3^/a^) 

= -2q;3-2(v^ + q;3V^) < 0. (20) 

From Eqs. (11), (18) and (19) and by using t{z) — f{z) + 2t:{z) and equating Eq. (11) 
with the condition that A must satisfy given by Eq. (12), we obtain the energy eigenvalue 
equation for the potential under the consideration 

n[{n-l)a2, + a-2-2a^\ - ^ {2n ^ l){^/ol^ + azy/oi^) 

+ q;7 + 2q;3Q;8 + 2y/a^ag = . (21) 

By using Eq. (9) 

p(^) = ^«io-i(i_ ^3^)^-10-1^ (22) 
and together with Eq. (8), we obtain 

^„(.) = Pi""-^' - -""-^1 - 2.3.) , (23) 

where ctio = cti + 2q;4 + 2y^ , ctn — a'z — 2ol^ + 2(^^09 + ol-^^Jol^^ and P^'^^{\ — 2asz) are 
Jacobi polynomials. By using Eq. (7), we obtain 

^jj(z) = (1 - a^zy''-^ , (24) 
and the total wave function become 
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where 0:12 = ^4 + ^/^s j <^i3 = ct5 ^ (^/o^g + ^^s^/f^s)- 

In some problems the situation appears where eta = 0. For this type of the problems, the 
solution given in Eq. (25) becomes as 



^>{z) = ^"i^e"""L^i°-^(aii^) , (26) 

and the energy spectrum is 

— 2a^n + {2n + 1)( ^Jo^^ — a^^/as) + n{n — 1)0:3 + ct? 

+ 2q;3Q;8 - 2^Jasag + 0:5 = . (27) 

when the limits become lima^^oPn ' {l — a^z) — L'^^'^~^ {auz) and ZimQ3_>o(l — 
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IV. BOUND STATES 

A. The Hypergeometric-Type Potential 

The Dirac equation has the exact pseudospin symmetry if E(r) = C = const. , so Eq. (3) 
becomes under that condition 



- ^^7^ -{m-E + C)MA(r)}G,,(r) = , (28) 

where «; = ^+ lforK>0, and k — —I lor k < Q . The hypergeometric-type potential is 
given by 



■ n _|_ n p-2ar \ 2 

F(r) = D[l - a coth(ar)]^ = ( \ _ ^!,,, ) , (29) 



where the real parameters D, cr, and a represent the potential [22], and Di = v-D (1 — cr) , 
and D2 = v^(l + (7). 



Eq. (28) can not be solved analytically for any k values because of k{k — l)/r^ term, so 
we use the approximation ~ 4Q;^e"^"''/(l — e~^"'')^ [26] to solve the equation for any 
spin-orbit quantum number k . 

By using this approximation to centrifugal-like term, setting A(r) to the potential given 
in Eq. (29), and inserting into Eq. (28), we obtain 



where fi = m — E + C , and e = m{m + C) + E{C — E) . At this point, it is worthwhile to 
note that Eq. (29) becomes for a = l{Di = 0) 



= ^2 1^13^ J ' (31) 
This form of the potential corresponds to the Manning- Rosen potential for A = [27] if we 



set y ^ a(Q; — 1) ^ D2 (here, k and a are the parameters in Ref. [27]), and | — >• 2q; . This 
means that we could also obtain the energy eigenvalue equation of the Dirac equation for 
the Manning-Rosen potential in the case of the exact spin symmetry, if we set the parameter 
Di = in the equations. 

By using the new variable z — e~^"^''(0 < z < 1) , we obtain from Eq. (30) 



d^Gnniz) . 1 - Z dGnniz 



dz^ z{l-z) dz [zil-z)]^ 

+ 2(3^[DiD2i^ - e - 2a\{K - l)]z + p^i^Dj - e)z^}Gn^{z) = . (32) 

By comparing Eq. (32) with Eq. (13), we get the parameter set 





1, 


0.2 — 


1, 


0:3 = 


1, 


a4 = 


0, 


as = 


1 

2' 


ae = 


Ci + i 


0:7 = 


-6, 


as = 


6, 


aio = 




an = 


-2 + 21 


(X12 = 






1 

2 



5 = 2(3^[DiD2fi - e - 2a'^K{K - 1 



-^-(^6-6 + 6 + 1 +v^). 
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The energy eigenvalue equation becomes 



(Jp^l^ia^K - 1) - //(L>i + L>2)2] + I + /3^e - fxDf j (2n + 1 + 2/3^e - fxDf ) 

[4a\(fi: - 1) - 2(e + = -n(n + 1) - 1/2 . (34) 

In this case, we use only the negative energy eigenvalues, because negative energy states 
exist in the pseudospin symmetry [28]. 

Now, let us give the corresponding Dirac spinors from Eq. (25) 



X 



(35) 

Finally, we briefly give the energy eigenvalue equation for the special case a — 1 , which 
gives the energy spectra of the Manning-Rosen potential with A — 



(j4:f3^[a^K{K-l)- Dij] + ^ + j (2n + 1 + 2/3^/^ ) 



+4/3^ a'^K{K - 1) 



+ 4((2n + l)' + l) = 0. 



Under the exact spin symmetry, i.e. A(r) — C — const. , Eq. (4) becomes 



(36) 



f (P k{k + 1) 



m + E-C)M^{r)}Fn.{r)^0, 



(37) 



where K = —{£+ 1) for k > and k = £ for k < . By setting S(r) to the potential given 
in Eq. (30), using the above approximation to the k{k + l)/r^ term, and using the new 
variable z — e~^"''(0 < 2; < 1), we obtain 



dz^ 



+ 



:{p\e'-n'Dl) 



1- Z dFnniz) 

z{l-z) dz [z{l-z)Y 
p^[2e' + 2DiD2ii' + ia\{K + 1)];^ + (5\e' - ii'Dl)z^}F„^{r) = 



(38) 
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where n' ^m + E-C^e' ^ m{C -m) + E{E -C), and /^^ = 1/4q;2. By comparing Eq. (39) 
with Eq. (13), we obtain the parameter set given in Eq. (34) where = — /3^(e' — A*'-D|), 
5 = -2/32[e' + DiD2i^' + 2a^K{K + 1)], and {3 = -/^^(e' - ^^'Dl). The energy eigenvalue 
equation of the hypergeometric potential for the exact spin symmetry is written from Eq. 
(21) as 



(Jf3^li'{D^ + A)^ + + + f3^fx'Dl - e' j (2n + 1 + 2f3^fj,'Dj - e' ) 

[Aa^K{K + 1) + 2(e' + At'L>iL>2)] = -n{n + 1) - 1/2 , (39) 

The last equation can give negative, and positive eigenvalues, but we choose only positive 
energy eigenvalues, because in the case of the exact spin symmetry appears only the positive 
energy eigenstates [28]. The corresponding Dirac spinors are obtained from Eq. (25), and 
given 



(40) 

Now, we briefly give the energy eigenvalue equation for the Manning-Rosen potential 
with A^O (a^l) 



(j4:PW{m + E-C) + k{k + 1) + ^ + (3^m{m - C) + E{C - E)^ 

X (2n + 1 + 2(3^Jm{m - C) + E{C - E) ) 

+ /^^ Ua^K{K + 1) + 2(m(C - m) + E{E - C))] + - {{2n + 1)^ + 1) = . (41) 



B. The Generalized Morse Potential 



Assuming that the potential A(r) = Vy{r) + Vs{r) is the generahzed Morse potential given 



by 



V{r) = Fie-^"'' - V^e 



(42) 
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and substituting Eq. (42) into Eq. (3), and taking E(r) — 'E — const., we have the following 
equation in the exact pseudospin symmetry ior k — {z — e""'') 



{ 



+ 



1 d 



dz^ s dz 



+ 



4/3'(At'-£;2 + £;(// + E)) 



+ 4/3V2(/x -E + Y.)z + A(5^Vi{E -fx- J:)z^] }Gnniz) = . (43) 
Comparing Eq. (43) with Eq. (13), we obtain the following parameter set 



«! = 1 , 
as = 0, 
0(5 = 0, 
Oi7 = -6, 



a2 = 0, ei = 4/3Vi(/i-E + S), 

as = 6 , ag = 6 



"12 = \/|3, Q;i3 = -v^- 

Prom Eq. (27), we obtain the energy eigenvalue equation for k = 



(44) 



E'-E{f, + E)-fx'^^^(2n+l-^^l,-E + E^ , (45) 

where = l/4«^. Wc should choose the negative energy sohition in Eq. (46) because the 
negative energy states exist only in the exact pseudospin limit. The corresponding lower 
spinor component can be obtained from Eq. (26) 



X (4/J^vi(^_£; + s) .) , (46) 

In the case of exact spin symmetry the potential A(r) = Vv{r) — Vs{r) is a constant, let 
say A(r) = A = const.. We set the potential E(r) as the Morse potential in Eq. (42). 
Substituting the potential into Eq. (4), and using the same variable z — e~"'', we obtain 



+ Ap^V2{pi + E-A)z + Ap^V,{-E-pi + A)z^]}F^,{z)^0, (47) 



11 



Comparing Eq. (47) with Eq. (13), we obtain the parameter set given in Eq. (44) where 
^1 = Ap^Vi{E + At - A), ^2 = ^P^V2{E + 11- A) and ^3 = 4/32(^,2 _ + A{E - /i)). Prom 
Eq. (27), we obtain the energy eigenvalue equation in the case of exact spin symmetry for 

K = 



/.^ - + A(E - /.) = ^ (2n + 1 - ^ VW^) ' • (48) 

where = 1/4q;^. We should choose the positive energy solution in Eq. (49) because the 
positive energy states exist only in the exact spin limit. The corresponding Dirac spinor can 
be written as 



X ^^f^^y^^^^E-A) z) . (49) 



C. The Poschl- Teller Potential 



By taking the potential A(r) = Vy{r) + Vs{r) is the Poschl- Teller potential [23] given by 



-2ar 

= -4H^YT^z^ , (50) 

and substituting Eq. (50) into Eq. (3), taking into account E(r) — "E — const., we have the 
following equation in the exact pseudospin symmetry ior k — {z — —e~^°''^) 



dz'^ 



+ 



z{l — z) dz 
+ {ii + E)z'^]}G^,{z) = Q. 



H + E -[2H + 2E + AVo\z 



(51) 



PoUowing the same procedure, we obtain the parameter set 
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«! = 


1 


^1 = 




a2 = 


-1 

1, 


6 = 


2jJ''{-fi - h + E)[n + E + 2Vo\ , 


as = 


-1 

1, 






a4 — 


0, 


as = 


1 

2' 


0(6 = 


6 + 


q;7 = 




as = 




ag = 


= ^1 — ^2 + ^3 + 4, 


ttlO = 




an = 


= 2 + 2(y/6-6 + e3 + i 


"12 = 




"13 = 


= -|-(V^i-6 + e3 + i +\/!^ 



(52) 



and the energy eigenvalue equation of Poschl- Teller potential under the exact pseudospin 
symmetry for k = from Eq. (21) 



- -Y,{ii^E) = ^((2n+l)Q;+^4yo(/i-£^ + S) + Q;2y _ (53) 

The last energy eigenvalue equation has a quadratic form in terms of energy E. We take 
the negative energy vales in the exact pseudospin limit. The corresponding Dirac spinor can 
be written in terms of Jacobi polynomials, i.e., P^'^'^^x), 



In the case of exact spin symmetry, we set the potential E(r) as Poschl-Teller potential 
given in Eq. (50), A(r) = A = const., and by using the coordinate transformation z — 
_^-2ar ^ we obtain from Eq. (4) 



dz^ z{l — z) dz 

+ fn -E)- 2P'{pi -E- 2V,)z 

+ (5\ii-E)z']F^,{z) = Q, (55) 

which gives the parameter set given in Eq. (52) where ^1 = P'^{E — iJ,){A — fj, — E), ^2 = 
2f3^{A - n - E){E - n + 2Vo) and ^3 = f3'^{E - i^){A - ^ - E). The energy eigenvalue 
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equation of the Poschl-Teller potential under the exact spin symmetry ior k, — from Eq. 
(21) is obtained 

+ 1^'' - A{i^ - E) = ^(^{2n+l)a+^Wo{i^ + E-A)+a^y. (56) 
Finally, we obtain the corresponding Dirac spinor from Eq. (25) 



V. CONCLUSION 



We have studied the energy eigenvalues and the corresponding eigenfunctions of the Dirac 
equation with the hypergeometric potential, the Morse potential, and the Poschl-Teller po- 
tential in the case of pseudospin, and spin symmetry. We have used the parametric general- 
ization of the NU method to obtain the results. The energy eigenvalues of all potentials are 
real and the wave functions are written in terms of the Laguerre (Jacobi) polynomials. We 
have also investigated the special case a = 1 in the case of the hypergeometric potential, 
which corresponds to the case of the Manning-Rosen potential with /I = . So we have 
obtained the energy eigenvalue equation of the Manning-Rosen potential in the pseudospin, 
and spin symmetry case, respectively. 
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